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REGULARITY OF C 1 SURFACES WITH PRESCRIBED MEAN CURVATURE 
IN THREE-DIMENSIONAL CONTACT SUB-RIEMANNIAN MANIFOLDS 

MATTEO G AT. T.T AND MANUEL RITORE 


Abstract. In this paper we consider surfaces of class C 1 with continuous prescribed mean curvature 
in a three-dimensional contact sub-Riemannian manifold and prove that their characteristic curves 
are of class C 2 . This regularity result also holds for critical points of the sub-Riemannian perimeter 
under a volume constraint. All results are valid in the first Heisenberg group H 1 . 


1. Introduction 


Recently Cheng, Hwang and Yang [9], [10], have considered the functional 
(1.1) &(u) = f |Vii + F|+ f fu, 

on a domain fi c R 2 ", where F is a vector field and / e L°°(fi). In case F(x,y) = (— y, x), the 
integral f n |Vu + F\ is the sub-Riemannian area of the horizontal graph of the function u in the 
Heisenberg group H n . Among several interesting results, they proved in [10, Thm. A] that, in 
case n = 1, u e C 1 (fl) is an stationary point of & and / e C°(n), the integral curves of the 
vector field ((Vu + Fj/lVu + Fl) 1 , defined in the set |Vu + F| ^ 0, are of class C 2 . The geometric 
meaning of their result is that the projection of the characteristic curves of the graph of u are of 
class C 2 . A stationary point u of JP satisfies weakly the prescribed mean curvature equation 


( 1 . 2 ) 


div 


Vu + F A = 

|Vu + F\J ~ 


Theorem A in [10] is well-known for C 2 minimizers and generalizes a previous result by Pauls 
[21, Lemma 3.3] for H-minimal surfaces with components of the horizontal Gauss map in the 
class W 1,1 . For lipschitz continuous vanishing viscosity minimal graphs, it was proven by Ca- 
pogna, Citti and Manfredini [4, Cor. 1.6]. 

In order to extend this result to arbitrary surfaces, it is natural to replace & by the sub- 
Riemannian prescribed mean curvature functional 


(1.3) /(£,B) = P(£,B) + f /, 

JEnB 

where E is a set of locally finite sub-Riemannian perimeter in fi, P(£,B) is the relative sub- 
Riemannian perimeter of E in a bounded open set B c fi, and / e L°°(f2). If F c HP is the 
subgraph of a function t = u(x,y) in the Heisenberg group H", then ^(£) coincides with (1.1) 
taking £(x,y) = (— y, x). The notion of sub-Riemannian perimeter used in sub-Riemannian 
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geometry was first introduced by Capogna, Danielli and Garofalo [5] for Carnot-Caratheodory 
spaces. General properties and existence of sets with minimum perimeter were proved later by 
Garofalo and Nhieu [17]. A rather complete theory of finite perimeter sets in the Heisenberg 
group IF following De Giorgi’s original arguments was developed by Franchi, Serapioni and 
Serra-Cassano [11], and later extended to step 2 Carnot groups [12] by the same authors. The 
recent monograph [6] provides a quite complete survey on recent progress on the subject. 

We have defined the prescribed mean curvature functional following Massari [19], who consid¬ 
ered minimizers of J? for the Euclidean relative perimeter. He obtained existence and regularity 
results for this problem and observed that, in case E is the subgraph of a Lipschitz function u 
defined on an open bounded set D c R n_1 , the function u satisfies weakly the prescribed mean 
curvature equation 

div f --= : (x)=/U,u(») 

V^/l + l Vu\ 2 J 

for x e D. In case dE n El is a hypersurface of class C 2 then the mean curvature of dE at a point 
p € dE equals g(p). See also Maggi [18, pp. 139-140]. 

The aim of this paper is to extend Cheng, Hwang and Yang’s regularity result for characteristic 
curves [10, Thm. A] from C 1 horizontal graphs satisfying weakly the mean curvature equation in 
the first Heisenberg group HI 1 to surfaces of class C 1 with prescribed mean curvature in arbitrary 
three-dimensional contact sub-Riemannian manifolds. 

In this setting, the Euclidean perimeter is replaced by the sub-Riemannian one and the in¬ 
tegral of the function / is computed using Popp’s measure [20, § 10.6], [2]. The minimizing 
condition will be replaced by a stationary one. Our ambient space will be a three-dimensional 
contact manifold with a sub-Riemannian metric defined on its horizontal distribution. In partic¬ 
ular, no assumption on the existence of a pseudo-hermitian structure is made. We shall prove in 
Theorem 4.1 

Let E c El be a set with C 1 boundary and prescribed mean curvature / e C°(f2) 
in a domain El c M of a three-dimensional contact sub-Riemannian manifold. 

Then characteristic curves in dE are of class C 2 . 

We remark that [10, Thm. A] states that the projection of characteristic curves to the plane t = 0 
is of class C 2 , but together with [8, (2.22)] this implies that the characteristic curves themselves 
are C 2 . We thank J.-H. Cheng for pointing out this fact. 

While the proof of [10, Thm. A] was based on the integral formula [10, (2.3)], see also (3.7) 
in [16, Remark 3.4], the proof of Theorem 4.1 is purely variational and follows by localizing the 
first variation of perimeter along a characteristic curve. A much weaker version of Theorem 4.1 
was given in [16, Thm. 3.5], where it was proven that the regular part of an area-stationary sur¬ 
face of class C 1 in the sub-Riemannian Heisenberg group H 1 is foliated by horizontal geodesics. 
Theorem 4.1 provides a new result even for the case of the first Heisenberg group H 1 . 

The regularity of characteristic curves proven in Theorem 4.1 allows us to define in Section 5 a 
mean curvature function H in the regular part of dE, that coincides with /. As a consequence of 
the definition of the mean curvature, we shall prove in Proposition 5.3 that characteristic curves 
are of class C k+2 in case / is of class C k when restricted to a characteristic direction. This holds, 
e.g., when / e C k (El) of / e Cg(f2), the space of functions with continuous horizontal derivatives 
of order k, k ^ 1. This class contains C 1 (Q) when k ^ 2. Critical points of the perimeter, even¬ 
tually under a volume constraint, and C 1 boundary, have constant prescribed mean curvature as 
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shown in Section 3. Hence Theorem 4.1 applies to these sets and implies that the regular parts 
of their boundaries are foliated by C 00 characteristic curves, see Proposition 5.4. 

We have organized this paper into several sections. In the second one we provide the necessary 
background on contact sub-Riemannian manifolds and sets of finite perimeter, and we recall the 
first variation formula for C 1 surfaces following [13]. In Section 3 we introduce the definition of 
set of locally finite perimeter with prescribed mean curvature and prove that a set with C 1 bound¬ 
ary and area-stationary under a volume constraint has constant prescribed mean curvature. The 
main result, Theorem 4.1, is proven in Section 4. The consequences on the mean curvature and 
higher regularity for characteristic curves will appear in Section 5. 

2. Preliminaries 

2.1. Contact sub-Riemannian manifolds. In this paper we shall consider a 3-dimensional C°° 
manifold M with contact form to and a sub-Riemannian metric g^ defined on its horizontal dis¬ 
tribution 34? := ker(to). By definition, dto \ M is non-degenerate. We shall refer to (M,to, g^) as 
a 3-dimensional contact sub-Riemannian manifold. It is well-known that to A dto is an orientation 
form in M. Since 

dto(X, 7) = X(co(Y')) - 7(to(X)) - coax, Y ]), 

the horizontal distribution 34? is completely non-integrable. The Reeb vector field T in M is the 
only one satisfying 

(2.1) to(T) = 1, S? T co = 0, 

where if is the Lie derivative in M. 

A canonical contact structure in Euclidean 3-space R 3 with coordinates (x,y, t ) is given by 
the contact one-form ro 0 := dt + xdy — ydx. The associated contact manifold is the Heisenberg 
group H 1 . Darboux’s Theorem [3, Thm. 3.1] (see also [15]) implies that, given a point peM, 
there exists an open neighborhood U of p and a diffeomorphism cp p from U into an open set of R 3 
satisfying 4>*co 0 = co. Such a local chart will be called a Darboux chart. Composing the map f p 
with a contact transformation of H 1 also provides a Darboux chart. This implies we can prescribe 
the image of a point pel/ and the image of a horizontal direction in T p M. 

The metric g x can be extended to a Riemannian metric g on M by requiring T to be a unit 
vector orthogonal to 34?. The Levi-Civita connection associated to g will be denoted by D. The 
integral curves of the Reeb vector field T are geodesics of the metric g. This property can be easily 
checked since condition if r a> = 0 in (2.1) implies a>([T,X]) = 0 for anyX e 34?. Hence, for any 
horizontal vector field X, we have 

g(X,D T T) = -g{D T X, T) = -g(D x T, T) = 0. 

We trivially have g(T,D T T) = 0, and so we get D t T = 0, as claimed. 

The Riemannian volume element in (M,g) will be denoted by dM. It coincides with Popp’s 
measure [20, § 10.6], [2]. The volume of a set E c M with respect to the Riemannian metric g 
will be denoted by |£|. 

2.2. Torsion and the sub-Riemannian connection. The following is taken from [13, § 3.1.2]. 
In a contact sub-Riemannian manifold, we can decompose the endomorphism X e TM —» D X T 
into its antisymmetric and symmetric parts, which we will denoted by J and t, respectively, 

2 g(J(X),Y) = g(D x T,Y ) - g(D Y T,Xf 
2g(z(XlY) = g(D x T,Y) + g(D Y T,Xf 


( 2 . 2 ) 
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Observe that J(X), t(X) e #€ for any vector field X, and that J(T ) = t(T) = 0. Also note that 

(2.3) 2 g(J(X),Y) = -g([X,Y],T), X,Y£3f. 

We will call t the ( contact ) sub-Riemannian torsion. We note that our J differs from the one 
defined in [14, (2.4)] by the constant g{\X, Y], T), but plays the same geometric role and can be 
easily generalized to higher dimensions, [13, § 3.1.2]. 

Now we define the ( contact ) sub-Riemannian connection V as the unique metric connection, 
[7, eq. (1.5.3)], with torsion tensor Tor(X, Y) = Y7 X Y — V Y X — [X, Y] given by 

(2.4) Tor(X, Y) := g(X, T) T(Y) - g(Y, T) z(X) + 2g(J(X), Y) T. 

From (2.4) and Koszul formula for the connection V it follows that T is a parallel vector field for 
the sub-Riemannian connection. In particular, their integral curves are geodesics for the connec¬ 
tion V. 

If X e Jt?, p e M, and X p / 0, then J(X p ) ^ 0: as dco\ x is non-degenerate, there exists Yeii" 
such that dco p (X p , Y p ) ^ 0. From (2.2) we have 2g(J(X p ), Y p ) = — g([X, Y] p , T p ), different from 
0 since a> p ([X, Y] p ) = -deo(X p , Y p ) / 0. 

The standard orientation of M is given by the 3-form co A dco. If X p is horizontal, then 
the basis {X p ,J(X p ),T p } is positively oriented. To check this, observe first that the sign of 
(to A dco)(X, J{X), T) equals the sign of dco{X,J{X)), and we have 

doo(X,J(X)) = -co([X, J(X)]) = -g([X, J(X)], T) = g(Tor(X, JQ0), T) = 2g(J(X), J(X)) > 0. 

2.3. Perimeter and C 1 surfaces. A set E c M has locally finite perimeter if, for any bounded 
open set B c M, we have 

P(E,B):= sup] / divNdM : 17 horizontal, supp(I7) C B, IlLilloo ^ 1 l <+oo. 

[Jeob J 

The quantity P{E,B ) is the relative perimeter of E in B. 

Assuming S = dE is a surface of class C 1 , the relative perimeter of E in a bounded open set 
B c M coincides with the sub-Riemannian area of S n B, given by 

(2.5) A(EnB)= / \N h \dT,. 

inns 

Here N is the Riemannian unit normal to E, N h is the horizontal projection of N to the horizontal 
distribution, and dE is the Riemannian area measure, all computed with respect the Riemannian 
metric g, see [6]. The quantity \N h \ vanishes in the singular set E 0 c S of points peE where 
the tangent space T p E coincides with the horizontal distribution J^ p . The horizontal unit normal 
at p e E \ E 0 is defined by (v h ) p := (JV h ) p /|(iV fl ) p |. At every point peE \ E 0 , the intersection 
Jtf’p n T p E is one-dimensional and generated by the characteristic vector field Z := J(v h )/|J(v h )|. 
The vector S p is defined for p e E \ E 0 by S p := g(_N p , T p )(v h ) p — |(iV h ) p | T p . The tangent space 
T p E, p e E \ E 0 , is generated by { Z p ,S p }. 

2.4. The first variation of the sub-Riemannian perimeter for C 1 surfaces. Given a set E with 
C 1 boundary, we can use the flow {(p s ] seR of a vector field U with compact support in B to pro¬ 
duce a variation of E DB. The Riemannian area formula gives the following expression of the 
sub-Riemannian area of E s := <^ S (E n B), 

A(E s ) = ^|N ft s | Jac(^)dE, 

where N s is a unit normal to E,. Fix p e E \ E 0 and the orthonormal basis {e l5 e 2 } = ]Z p ,S p ] in 
T p E. We consider extensions E lt E 2 of Z p , S p , respectively, along the integral curve of U passing 
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through p. The vector fields E-^s), E 2 (s) are invariant under the flow of U and generate the 
tangent plane to E s at the point <p s (p). The vector (E-l x E 2 )/\{E l x E 2 )| is normal to E s . Here x 
denotes the cross product with respect to a volume form 17 for the metric g inducing the same ori¬ 
entation as co/\dco, i.e. g(w,uxv) = r](w,u,v ). It is easy to check that KEj xE 2 )|(s) = Jac(ip s )(p), 
and that 


Hence 


and we get 


V(p,s) := (Ej x E 2 ) h (s) = (g(E l5 T) (T x E 2 ) - g(E 2 , T) (E : x T))(s). 
A(Ej = ^|V(p,s)|dS(p), 
g(VnV,^) 


d 

ds 


lv(s, p)l = 


5=0 


IV P 


Since {(v h ) p ,Z„, T p } is positively oriented, observe that V„ = |(A/ h ) p | (v h ) p . On the other hand, 


VuV = g& Up Ei,T p KT p x(.E 2 ) p )-g(y Ur E 2 ,T p X(E 1 ) p x T p ) 

-g-((^ 2 )p, 7 ’ P )CV c/p E 1 x Tp), 

and so 

j(V[i V, Vp) 

--= -g(V Up E 2 ,r p )+ |av h ) p |g(V Up E a x T p , (v h ) p ). 

Since 


g(V [/p E 2 , T p ) = g(V (£2)p U + Tor(f/ p , (E 2 ) p ), T p ) 
= S p (g(H, T)) + g(Tor([/ p ,Sp), T p ) 
= S p (g([/,T)) + 2g(J({/ p ),S p ), 


g(V Up E! x T p , (v h ) p ) = g((V (£i)p [/ +Tor([/ p , (EOp)) x T p ,(v h ) p ) 
= t?((v ft ) p , V (£i)p U + Tor(H p , (E^p), T p ) 

— +g(V Zp !7 + Tor( Up, Z p ),Z p ) 

= +g(v Zp f/,z p ) + g([/ p , r p )g(T(z p ),z p ), 

we conclude that the first variation of the sub-Riemannian perimeter is given by 


( 2 . 6 ) 


d 

ds 


A(S s ) 

s=0 



{ - S(g(l7, T)) - 2g(J(H),S) + |AI),| g(V z I7,Z) 


+ |N h |g(t/ J T)g(T(Z),Z)}dE. 


This formula was obtained in [14, Lemma 3.4]. 


3. Sets with prescribed mean curvature 

The reader is referred to [18, (12.32) and Remark 17.11] for background and references in 
the Euclidean case. Consider a domain ilcM, and a function / : Q —* R. We shall say that a set 
of locally finite perimeter E cfi has prescribed mean curvature f on El if, for any bounded open 
set B c El, E is a critical point of the functional 

P(E,B)- [ f, 

Jehb 


(3.1) 
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where P{E,B ) is the relative perimeter of E in B, and the integral on £ nB is computed with 
respect to the canonical Popp’s measure on M, see [20] and [2], The admissible variations for 
this problem are the flows induced by vector fields with compact support in B. 

If E = dE is a surface of class C 1 in fi, then E has prescribed mean curvature / if it is a critical 
point of the functional 

(3.2) A(SnB)-/ /, 

Jehb 

for any bounded open set B cU, 

If £ is a critical point of the relative perimeter P(£,B) in any bounded open set Ben, then E 
has zero or vanishing prescribed mean curvature. 

Assume now that E c is a set of locally finite perimeter with C 1 boundary E, and that E is a 
critical point of the perimeter under a volume constraint. This means (d/ds) s= 0 A(ip s (EnB)) = 0 for 
any flow associated to a vector field with compact support in H satisfying (d/ds) s= 0 |i^ s (BnB)| = 0. 
If the perimeter of E in Q is positive, then there exists a (horizontal) vector field U 0 with compact 
support in n so that / Bnn div U 0 dM > 0. By the Divergence Theorem, 


f g(U 0 ,N) dE / 0, 
J Enf2 


where N is the outer normal to E. Let {t/bLem t> e the flow associated to the vector field U 0 and 
define 


(3.3) 


H 0 := 


d/ ds| s=0 A(i/> s (E)) 


d/ds\ s=Q \xp s (E)\ 

Let B c 11 be a bounded open subset and W a vector field with compact support in B and associ¬ 
ated flow {ip s } seK . Choose Aetso that W — A U 0 satisfies 


d 

ds 


s= 0 


|^(£)|-A — 


IVb(£)l= / S(W — hU 0 ,N)dT, = 0. 
=o ■/£ 


Then the flow associated to W — A U 0 preserves the volume of E n (B u B 0 ), where B 0 c Q is a 
bounded open set containing supp(D 0 ). Let Q(t/) be the integral expression in (2.6). From our 
hypothesis and the linearity of (2.6), Q(W — AI/ 0 ) = 0. Hence Q(W) = AQ(I/ 0 ). From (3.3) we 
get 


Q(W) = AQ(l/ 0 ) = AH 0 — 
ds 


s =0 


|^(£)|=H 0 — 


l‘Ps(£)l, 

s =0 


and so E has (constant) prescribed mean curvature H 0 


4. Main result 


In this Section we shall prove our main result 


Theorem 4.1. Let M be a 3-dimensional contact sub-Riemannian manifold, H C M a domain, and 
E c n a set of prescribed mean curvature f e C°(f2) with C 1 boundary E. Then the characteristic 
curves in E are of class C 2 . 


Proof. Given any point peE \ E 0 , consider a Darboux chart (U p , such that <£ p (p) = 0. The 
metric g x can be described in this local chart by the matrix of smooth functions 


G = 


S n S 12 
S 21 S22 


g(.x,x ) g(x,y) 

g(Y,X) g(Y,Y) 
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After a Euclidean rotation around the t-axis, which is a contact transformation in H 1 [22, p. 640], 
we may assume there exists an open neighborhood BnEofpeE\ E 0 , where B c HI 1 is an open 
set containing p, so that B n E is the intrinsic graph G u of a C 1 function u : D —» R defined on a 
domain D in the vertical plane y = 0. We can also assume that E n B is the subgraph of u. The 
graph G u can be parameterized by the map f u : D —> R 3 defined by 

f u (x, t) := (x, u(x, t), t - xu(x, t)), (x, t) € D. 


The tangent plane to any point in G u is generated by the vectors 

j~ x -» (1 ,u x , -u - xu x )=X + u x Y - 2uT, 
j- t -» (0,u t , 1 — xu t ) = u t Y + T, 

and so the characteristic direction is given by Z = Z/\Z\, where 

Z = X + ( u x + 2uu t )Y. 

If y(s) = (x(s), t(s)) is a C 1 curve in D, then 

T(s) = (x(s), u(x(s), t(s)), t(s) - x(s)u(x(s), t(s))) c G u 

is also C 1 , and so 

r'(s) = x' (X + u x Y - 2uT) + t' ( u t Y + T) = x'X + (x'u x + t'u t )Y + (t' - 2ux')T. 

In particular, horizontal curves in G u satisfy the ordinary differential equation t' = 2 ux'. Since 
u e C 1 (D), we have uniqueness of characteristic curves through any given point in G u . 

A unit normal vector to E is given by JV/[JV|, where 

N = {X + u x Y - 2uT) x (u t Y + T). 


Here x is the cross product with respect to the Riemannian metric g and a given volume form 
Tj chosen so that 17 (X,Y, T) > 0. Hence g(w,u x v) = i 7 (w, u, v). If { 62 , 62 , 63 } is an orthonormal 
basis so that p^e 1 ,e 2 ,e 3 ) = 1 and A is the matrix whose columns are the coordinates of X, Y, T 
in the basis {e l5 e 2 ,e 3 }, then r](X, Y,T) = det(A). On the other hand, as 



r gn 

§12 

0 ^ 

A l A = 

§21 

§22 

0 


V 0 

0 

V 


we get det(A ) 2 = det(G). Since det(A) > 0 we obtain det(A) = det(G ) 1|/2 and so 

r)(X,7,T) = det(G) 1/2 . 


Let E 1 = X + u x Y — 2uT, E 2 = u t Y + T. We compute the scalar product of N = E 1 x E 2 with X, 
Y, T to obtain 


g(X,E 1 xE 2 ) 

g{Y,E 1 xE 2 ) 

g(T, E 1 x E 2 ) 





f 1 

1 

0^ 




,e 2 ) 

= det 

0 

u x 


V(X, Y,T) = 

(■ u x + 2uu t ) det(G ) 1 ^ 2 




Vo 

-2 u 







f 0 

1 

0^ 



ri(Y,E 1 

e 2 ) = 

= det 

1 

U x 


rj(x,y,r) = 

— det(G) 1/2 . 





-2 u 







f 0 

1 

0^ 



V( T , E i 

£2) = 

= det 

0 



r)(X, Y,T) = 

u t det(G) 1 ^ 2 . 




V 1 

-2 u 

u 
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Since g(Y, E 1 x E 2 ) < 0 and E n B is the subgraph of u, the vector field E 1 x E 2 points into the 
interior of E. If N = E 1 x E 2 = aX + fiY + jT, then 


whence 



§12 

0^ 

cA 


Fu x + 2 tap 

§21 

822 

0 

P 

= det(G) 1/2 

-1 

V 0 

0 

V 



V u t J 


y = det(G) 1/2 u t . 


Let us compute now the sub-Riemmanian area of the intrinsic graph G u . It is easy to check 
that dS = | E 1 x E 2 \dxdt, i.e., that Jac(/ U ) = \E 1 x E 2 \. Since \N h \ = x E 2 ) h \/\E l x E 2 \, then 
using (4.1) and the explicit expression of the inverse matrix G -1 we get 

\N h \ Jac(/J = Pi x£ 2 ) h |= ((a |8) G 

= (CS "22 °/JK + 2 uu t ) 2 + 2 (g 12 of u Xu x + 2uu t ) + (g n °/J) 1/2 . 



Finally, from (2.5) we obtain 

(4.2) A(GJ= [ (g 22 (u x +2uu c f+ 2g 12 (u x +2uu t ) +g n ) 1/2 dxdt, 

Jd 

where, by abuse of notation, we have written g i; - instead of the cumbersome notation (g i; - °/ u ). 

Now we consider variations of G u by graphs of the form s —» u + sv, where v e C“(D) and s 
is a real parameter close to 0. This variation is obtained by applying the flow associated to the 
vector field vY to the graph G u . The function v is obtained by extending v to be constant along 
the integral curves of the vector field Y, and multiplying by an appropriate function with compact 
support equal to 1 in a neighborhood of £. 

When F is a function of (x, y, t), we have 

f 8 F dF\ 

(F °fu+sv)(x,t)= 3- x— v(x,t) = y / CX;0 (F)v(x,t). 

s =o c)t J ux,t) 

So we get 

A(G u+iv )= / (JfjV + M (v x + 2uv t + 2vu t )) dxdt, 

s =0 D 

where the functions K 1 and M are given by 


d 

ds 


d 

ds 


K _ 1 Y(g 22 Xu x + 2uu t f + 2Y(g 12 )(tG + 2 uu t ) + Y(g n ) 
1 2 (g 22 (u x + 2uu £ ) 2 + 2g 12 (iz x + 2uiz t ) + g 11 ) 1/2 

and 

g 2 2(Ux+2uu t ) + g 12 

M =-;-. 

+ 2 uuj 2 + 2g 12 (u x + 2uu t ) + gn) 1/2 

Observe that the functions ffj and M are continuous. Since 


X + (u x + 2 uu t ) Y 

iS 22 ^ u x + 2 uu t ) 2 + 2 g 12 (u* + 2uu t ) + g n ) 1/2 ’ 


the function M coincides with g(Z, Y) °f u . A straightforward computation implies 


1 = |Z| 2 = det(G)- 1 (g 22 g(Z,Z) 2 - 2g 12 g(Z,X)g(Z, Y) + g n g(Z, Y) 2 ) 
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and so 


(4.3) 


g(Z,X) 


gi 2 g(Z, y) ± (det(G)(g 22 - g(Z, 7 ) 2 )) 1/2 
g22 


By Schwarz’s inequality g(Z, Y) 2 ^ g(Y,Y ) = g 2 2 - Inequality is strict since otherwise Y and Z 
would be collinear. Hence g(Z,X) has the same regularity as g(Z, 7) by (4.3). 

The subgraph of u can be parameterized by the map (x, t,s) —» (x,s, t — xs). The Jacobian of 
this map is easily seen to be equal to det(G). Hence 


d 

ds 


s=0 J subgraph G u _ 


f= f det(G)v dxdt. 


If E has prescribed mean curvature /, this implies 


(4.4) / (ffv + M(v x + 2uv t + 2vu t ) dxdt = 0, 

Jd 

for any v e C“(D), where the continuous function K is given by K = K 1 — f det(G). By Remark 
4.3 below, (4.4) also holds for any v e C°(D) for which v x + 2 uv t exists and it is continuous. 

Now we proceed as in the proof of Theorem 3.5 in [16]. Assume the point peG„ corre¬ 
sponds to the point (a, b ) in the xt-plane. The curve s —»(s, t(s)) is (a reparameterization of the 
projection of) a characteristic curve if and only if the function t(s) satisfies the ordinary differ¬ 
ential equation t'(s) = u(s, t(s)). For e small enough, we consider the solution t e of equation 
t'Js) = 2 u(s, t e (s)) with initial condition tja) = b + s, and define y e (s) := (s, t e (s)), with y = y 0 . 
We may assume that, for small enough e, the functions t e are defined in the interval [a — r,a + r] 
for some r > 0. The function dt e /de satisfies 


f dt e \ f dt e \ dt e 

(4.5) ( 77J M = 2u t (s,as))f jfjis), -gfW = !■ 

where ' is the derivative with respect to the parameter s. 

We consider the parameterization 


F^,e) :=(£,t,(O) = (s ,0 


near the characteristic curve through (a, b). The jacobian of this parameterization is given by 


det 


dtjde 


3t e 

de ’ 


which is positive because of the choice of initial condition for t e and the fact that the curves y f; (s) 
foliate a neighborhood of (a, b). Any function if can be considered as a function of the variables 
(£, e) by making ip(E,, e) := <p(£, t e (0)- Changing variables, and assuming the support of if is 
contained in a sufficiently small neighborhood of (a, b), we can express the integral (4.2) as 


<- a+r f (dif 

Kif + M l— + 2ifu t 


d t e 

-r d ^ 

de 


de. 


where I is a small interval containing 0. Instead of if, we can consider the function ifh/(t e+h — tj, 
where h is a sufficiently small real parameter. We get that 


d r hif h dif h _ u(J e + h) — u{E,, e) h 

d ? V1 e+h Ifi J d ^ ^e+h t e te+h t e te+h t e 


d(p/dE, 2 lfU t 
dtjde dtjde’ 


tends to 
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when h —> 0. So using that G u is area-stationary we have that 

JaL - T^T — )))j 7 d A ds 


vanishes. Furthermore, letting h —> 0 we conclude 


HL 


a+r f d(p 

K(p + M 




dE,\cLe = 0 . 


Let now 17 : R —» R be a positive function with compact support in the interval I and con¬ 
sider the family rj p (x) := p _ 1 rj(x/p). Inserting a test function of the form R p (.e)ip(f), where 
xp is a C 00 function with compact support in (a — r, a + r), making p —» 0, and using that G u is 
area-stationary we obtain 


f 


(W, O WO + M( 0, d^ = 0 


for any xp e C“((a — r, a + r)). By Lemma 4.2, the function M(0, which is the restriction of 
g(Z,Y ) to the characteristic curve, is a C 1 function on the curve. By equation (4.3), the restric¬ 
tion of g(Z,X) to the characteristic curve is also C 1 . This proves that horizontal curves are of 
class C 2 . □ 

Lemma 4.2. Let I c R be an open interval, k, m e C°(I), and K e C^J) be a primitive of k. 
Assume 


(4.6) 


J kxp + mxp' = 0 , 


for any xp e C“(/). Then the function —K + mis constant on I. In particular, m e C 1 (7). 
Proof. Since (Kxpf = kxp + Kxp', integrating by parts we see that (4.6) is equivalent to 

( - K + m) xp' = 0, 


/< 


for any xp e C“(7). This implies that —K + m is a constant function on I. 


□ 


Remark 4.3. Let us check that (4.4) holds for any w e C°(D) such that w x + 2uw t exists and is 
continuous. Let us consider a sequence w ; - e C“(D), where = pj *w, and p ; denote the stan¬ 
dard modifiers. We have that Wj converges to w and that (\Vj) x +2u(w J ) t converges to w x + 2uw t 
uniformly on compact subsets of D, for j —> 00. We conclude 

0 = lim / (.K(w;) + M((w.-)*+2u(w,-) t + 2w,u t ) dxdt = / (Kw + M(w x + 2uw t + 2wu t ) dxdt, 
.H 00 Jd Jd 

thus proving the claim. 

Remark 4.4. In case M is the Heisenberg group H 1 , G is the identity matrix and the expression 
for the sub-Riemannian area of the graph G u given in (4.2) reads 


A(G U ) = / ((i i x + 2uuf) 2 + 1 ) 1/2 dxd t, 


a well-known formula obtained in [ 1 ]. 
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5. The mean curvature for C 1 surfaces 


Given a surface E c M of class C 1 such that the vector fields Z and v h are of class C 1 along 
the characteristic curves in E \ E 0 , we define the mean curvature of E at p e E \ E 0 by 

(5.1) (div£(v h ))(p) := -g(V z v h ,Z)(p). 

This is the standard definition of mean curvature for C 2 surfaces, see e.g. [14, (3.8)] and the 
references there. The mean curvature is usually denoted by H. The mean curvature depends on 
the choice of v h . In case E is the boundary of a set E, we shall always choose N as the inner 
normal and v h = N h /\N h \. 

Using the regularity Theorem 4.1 we get 


Proposition 5.1. Let El cz M be a domain and E c El a set of prescribed mean curvature f e C°(fi) 
with C 1 boundary E with H e Lj oc (T,). Then the first variation of the functional (3.2) induced by a 
vector field U e Cq(£1) is given by 

(5.2) / H g(U,N)dT, — [ fg(U,N)dX. 

Jn Je 


Proof. We first observe that formula 


d 

ds 


M. E s ) 

s =0 


Hg(U,N)d E 


can be proved as in [14, Prop. 6.3], see also [14, Remark 6.4]. On the other hand, it is well- 
known that 


d 

ds 


s =0 


4.(n) 


/ = 


- J f g(U,N)dV, 


see e.g. [18, 17.8]. 


□ 


Then we have that the mean curvature H defined in (5.1) coincides with the prescribed mean 
curvature / 

Corollary 5.2. Let E C El be a set of prescribed mean curvature f e C°(f2) with C 1 boundary E in 
a domain El c M. Assume H <E L ; 1 oc (E). Then H[p) = f(p)for any p e E \ E 0 . 

Finally we can improve the regularity of Theorem 4.1 assuming the mean curvature function is 
more regular. This result is specially useful if we assume that higher order horizontal derivatives 
of the function / exist and are continuous. 

Proposition 5.3. Let E c El be a set of prescribed mean curvature f e C°(f2) with C 1 boundary E 
in a domain El c M. Assume that f is also C k in the Z-direction, k ^ 1. Then the characteristic 
curves of E are of class C k+2 inT,\ E 0 . 


Proof. Since we have V z Z = Hv h , we can write 

V z (V z Z) = Z[H)v h + HV z v h = Z{H)v h - H 2 Z. 

Iterating the procedure we obtain the statement. □ 

In particular, this result holds when / e Cg(fi), i.e., when / has horizontal derivatives of class 
(k — 1), see [11]. 

An important particular case is that of critical points of perimeter, possibly under a volume 
constraint. Assuming C 1 regularity of the boundary, these sets are known to have constant pre¬ 
scribed mean curvature from the discussion in Section 3. From Proposition 5.3, we immediately 
obtain 
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Proposition 5.4. Let E d LI be either a critical point of the sub-Riemannian perimeter or a critical 

point of the sub-Riemannian perimeter under a volume constraint. If E has C 1 boundary, then the 

regular part of dE is foliated by C°° characteristic curves. 
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